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Abstract 



■ Using BlackHat in conjunction with SHERPA, we have computed next-to-leading order QCD 



predictions for a variety of distributions in Z, 7* -|- 1, 2, 3-jet production at the Tevatron, where the 
Z boson or off-shell photon decays into an electron-positron pair. We find good agreement between 
the NLO results for jet pT distributions and measurements by CDF and DO. We also present jet- 
production ratios, or probabilities of finding one additional jet. As a function of vector-boson px, 
the ratios have distinctive features which we describe in terms of a simple model capturing leading 
logarithms and phase-space and parton-distribution-function suppression. 

PACS numbers: 12.38.-t, 12.38.Bx, 13.87.-a, 14.70.-e 
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I. INTRODUCTION 



The Large Hadron Collider (LHC) recently passed the milestone of first collisions. The 
start of the LHC era in particle physics opens new opportunities to confront data with 
theoretical predictions for Standard-Model scattering processes, at scales well beyond those 
probed in previous colliders. This confrontation will be a key tool in the search for new 
physics beyond the Standard Model. Where new physics produces sharp peaks, Standard 
Model backgrounds can be understood without much theoretical input. For many searches, 
however, the signals do not stand out so clearly, but are excesses in broader distributions of 
jets accompanying missing energy and charged leptons or photons. Such searches require a 
much finer theoretical understanding of the QCD backgrounds. 

An important class of backgrounds is the production of multiple jets in association with 
a Z boson. If the Z boson decays into neutrinos, this process forms an irreducible back- 
ground to LHC searches for new physics, such as supersymmetry, that are based on missing 
transverse energy and jets, as discussed in, e.g. ref. These processes can be cahbrated 
experimentally using events in which the Z boson decays into a pair of charged leptons, ei- 
ther an electron-positron pair or a di-muon pair. The latter samples are quite clean, and the 
QCD dynamics is of course identical to the Z ^ vv mode. The one experimental drawback 
is the small branching ratio for Z — )■ l^l^ . Nevertheless, there are already results from the 
Tevatron on Z production in association with up to three jets j2"5|, along with the prospect 
of new analyses using larger data sets in the near future. Therefore in this paper we focus 
on the production of Z + 1,2,3 jets at the Tevatron, which can also serve as a benchmark 
for future LHC studies. 

The first step toward theoretical control of QCD backgrounds at hadron colliders is the 
evaluation of the cross section at leading order (LO) in the strong coupling, as- Several 
computer codes are available for LO predictions. These codes typically use matching 
(or merging) procedures [ol, [lo| to incorporate higher- multiplicity leading-order matrix ele- 
ments into programs that shower and hadronize partons U lsl. Although such programs 
provide a hadron-level description which has great utility, the LO approximation suffers from 
large factorization- and renormalization-scale dependence, which grows with increasing jet 
multiplicity. This dependence is already up to a factor of two in the processes we shall study; 
accordingly, LO results do not generally provide a quantitatively reliable prediction. The 
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problems go beyond that of normalization of cross sections: shapes of distributions may or 
may not be modeled correctly at (matched) LO, and the results at this order can depend 
strongly on the functional form and value chosen for the scale. In order to resolve these prob- 
lems, and provide quantitatively reliable predictions, one must evaluate the next-to-leading 

y more 



r 



order (NLO) corrections to processes of interest. Such computations are technical 
challenging, but generically yield results with a greatly reduced scale dependence^ 
as well as displaying better agreement with measurements (see e.g. refs. 2|, y, [ 

More generally, sufficiently accurate QCD predictions can provide important theoretical 
input into experimentally-driven determinations of backgrounds, by allowing measurements 
in one process to be converted into a prediction for another, where theory is used only for the 
ratio of the two processes. For example, complete NLO predictions for W+ 3-jet production, 
followed by — )■ lu, are already available at parton-level 18|, and the present paper 



describes analogous predictions for Z+ 3-jet production, followed by Z — ^ These results 
allow the ratio oi Z to W production in association with up to three jets to be computed 
at NLO. Parton-level results do neglect non-perturbative effects, such as hadronization and 
the underlying event, which can contribute to both Z and W production. However, as 
long as the experimental cuts on the jets are the same, the non-perturbative corrections 
should largely cancel in the Z to W ratio. Therefore, measurements of the (more copious) 
production of W bosons in the presence of multiple jets can be extrapolated to the case of 
Z bosons 19|, using precise theoretical values for the ratio of Z to W events with similar 
kinematics, and as a function of the kinematics. Because leptonically decaying Z bosons, 
although rarer, are cleaner experimentally than W bosons, it has also been suggested to 
reverse t he p rocedure and use Z{^ l^l") boson samples to calibrate Iv) and Z{-^ vu) 

samples {20]. However, this procedure requires more input from theory: in order to have 
adequate statistics in the Z channel, less energetic kinematical configurations with 

larger cross sections are measured, and then extrapolated to more energetic configurations 
with smaller cross sections. 

In the past, the bottleneck in computing NLO QCD corrections to processes with large 
numbers of jets has been the evaluation of one-loop amplitudes involving six or more par- 
tons |l4]. On-shell methods 15, 21-30] have successfully resolved this bottleneck, by avoiding 



gauge-noninvariant intermediate steps and reducing the problem to much smaller elements 
analogous to tree-level amplitudes. In this paper we evaluate the required one-loop am- 
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32| . which implements on-shell 



phtudes with the BlackHat program hbrary llTI . Il8l |31 
methods numerically. A number of programs based on on-shell techniques have been con- 
structed by other groups [33 



371 . Approaches based on Feynman diagrams have also led 



to new results with six external partons, in particular the NLO cross section for producing 



four heavy quarks^ hadron colliders 
on-shell methods 



38 



. The ttbb case has also been computed via 



37j . and recently the final state tt + 2 jets has been computed at NLO 
in a similar way 

processes involving many external gluons, which often dominate multi-jet final states. 



NLO parton-level cross sections for the production o 
one or two jets have long been available in the MCFM 



a.W OT Z boson in association with 



amplitudes from ref. 



4l| code, which utilizes the one-loop 



22| for the two- jet case. More recently, complete NLO results have been 



obtained for W + 3-jet production [iSS] using BlackHat in conjunction with the SHERPA 
package [l3|. (Different leading-color approximations and "adjustment procedures" have 
also been applied to W + 3 jets at NLO 17, SgJ.) In this work, we use the same basic 
calculational setup as in ref. [l8| to compute Z + 3-jet production. The real emission. 



dipole subtraction 42| and integration over phase space is handled by AMEGIC-I — h 



which is part of the SHERPA package 
subtraction methods 



42 



13| . (Other automated implementations of infrared 



44j | have been described elsewhere [45|].) SHERPA is also used to 



perform the Monte Carlo integration over phase space for all contributions. One important 
improvement in the present study with respect to ref. [3] is to increase the efficienc y o f the 



phase-space integrator, making use of QCD antenna structures along the lines of refs. |46|,|47[. 

In this article we present results for Z, 7* + 3-jet production at the Tevatron to NLO 
in QCD, at parton level, and with the vector boson decaying to a lepton-anti-lepton pair. 
We include off-shell photon exchange, and 7* — Z interference, because the production of a 
charged-lepton pair by an off-shell photon is indistinguishable from the leptonic decay of a Z 
boson. Preliminary versions of some of the results presented here may be found in ref. 48| . 
where slightly different jet cuts were applied, along with a leading-color approximation. 

Here we present total production cross sections, with jet and lepton cuts appropriate to 
existing CDF and DO analyses [2, 4|, as well as a variety of distributions. We also study how 
Z, 7* -|- 3-jet production at the Tevatron depends on a common choice of renormalization 
and factorization scale /i. As mentioned above, LO results are generically rather sensitive 
to the scale choice. This sensitivity usually is greatly reduced at NLO. As an example, in 
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7* + 3-jet production, varying the scale by a factor of two from our default central value 
causes nearly a 60% deviation from the central value. In contrast, at NLO this deviation 
drops to only 15-22%. 

Total cross sections with standard experimental jet cuts are dominated by the production 
of jets with low transverse momentum, < My, where the vector boson V is a W or Z. 
Accordingly, scale choices such as the mass of the vector boson /i = My are reasonable for 
these quantities. For the study of differential distributions, however, it is better to choose 
a dynamical scale, event by event, in order to have reasonable scales for each bin 49] . This 
helps reduce the change in predicted shapes from LO to NLO, and can improve the NLO 
prediction some too. However, care is required in choosing the functional form of such 
scales. Greater care is required at the LHC than at the Tevatron, because the much larger 
dynamical range can ruin seemingly-reasonable scale choices, such as the commonly-used 
vector-boson transverse energy, /i = = J My + (p^)^ (see e.g. refs. 



As noted in ref. 



17 



50|). 



18| . a scale choice of /i = leads to negative cross sections in tails of some 



NLO distributions, because typical energy scales in the process are much larger than Ej, . In 
addition, as noted independently 5l|, the choice fi = E^ leads to undesirably large shape 
changes between LO and NLO. As we discuss here, for Z -\- 3-jet production at the Tevatron, 
the scale fi = E^ is unsatisfactory, at least at LO: choosing it results in large shape changes 



in distributions between LO and NLO. The tota^ 
fraction thereof), adopted in our previous study 



partonic transverse energy (or a fixed 



18|, is a satisfactory choice. Other choices. 



such as the combined invariant mass of the jets 5]J, should also be satisfactory. In any case, 
the lesson is clear: the vector-boson transverse energy is not satisfactory and should not be 
used, especially for processes at the LHC. 

At the Tevatron, both CDF and DO have measured {2, 4| jet-p^ distributions for 2', 7* + 
1, 2-jet production in the channel Z — )■ e~^e~. The DO measurements are not absolute cross 
sections, but are normalized to the inclusive Z, 7* -|- 0-jet cross section for the same set 
of lepton^ cuts. CDF and DO have each compared their data with NLO predictions from 
MCFM [41], taking into account estimates of non-perturbative corrections. For comparison 
we present our own NLO analysis of these processes. We then turn to the more involved case 
of 2', 7* -|- 3-jet production, which we compare against the DO experimental measurement. 
A difference between our NLO study and the experimental measurements is in our use of 
infrared-safe jet algorithms (as reviewed in ref. Our default choice here is the SISCONE 
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jet algorithm 55|], although we present some results using the anti-fcy algorithm j56| as well. 



(The kx algorithm 



gives parton-level results very similar to those for the anti-fcy 

le CDF and DO analyses, which 
59l-l6l|. are generically infrared 



one, so we do not present them.) The algorithms used in t 
are in the "midpoint" class of iterative cone algorithms 
unsafe, and cannot be used in an NLO computation of + 3-jet production. Although a 
midpoint algorithm will yield finite results at NLO for V + 2-jet production, it suffers from 
uncontrolled non-perturbative corrections that are in principle of the same order as the NLO 
correction 54l |. 



In comparing theory and experiment, the differing jet algorithms do introduce an addi- 
tional source of uncertainty. Nevertheless, based on our study of Z, 7* + 1, 2-jet production 
we expect the NLO results to match experiment reasonably well. There have also been two 
studies coniparing inclusive jet cross sections for midpoint algorithms with those for SIS- 



CONE 



55 



62], using Pythia [11], which find that the major non-perturbative differences 
between the algorithms are at the level of the underlying event. Although the kinematics 
of inclusive-jet production differs somewhat from that of a vector boson plus multiple jets, 
these results suggest that hadron-level data collected using SISCONE would differ from that 
for midpoint algorithms primarily by the underlying event correction (at least for larger 
cone sizes). That is, if the two measurements were corrected back to parton level, one would 
expect them to have only percent-level differences. (At the LHC, both ATLAS and CMS 
have adopted infrared-safe jet algorithms, removing this important source of uncertainty.) 

We shall present Z, 7* -|- 1, 2, 3-jet production cross sections for the CDF and DO cuts, as 
well as jet pr distributions. For Z, 7* -|- 3-jet production, we show the pr distributions for 
all three jets, ordered in p^- In addition, we discuss ratios of cross sections. Experimental 
and theoretical systematic uncertainties cancel to some degree in such ratios. Ratios of 
similar processes — for example, the ratio of -|- 3-jet to Z + 3-jet production — are 
thus attractive candidates for confronting experimental data with theoretical predictions. 
We shall not study this kind of ratio in the present paper, but another kind, that of a 
production process to the same process with one additional jet, sometimes known as the 



"jet-production ratio" 



63 



6J] . (This ratio is also called the "Berends" or "staircase" ratio. 



In particular, we study the ratio of Z, 7* -|- n-jet to Z, •y* + {n — l)-jet production up to 
= 4 at LO and = 3 at NLO. Its evaluation for different values of n can also test the 
lore, based on Tevatron studies, that the ratio is roughly independent of n. We find that for 
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total cross sections through n = 3, with the experimental cuts used by CDF, this scaling is 
valid to about 30%, and for the DO cuts, it is valid to about 15%. 

One may also ask whether not just total cross sections, but also differential distributions, 
can be predicted for 2', 7* + n-jet production (at least approximately) by scaling results 
for Z,'~f* + (n — l)-jet production. We show explicitly for n < 4, using the example of 
the vector-boson pr distributions, that shapes of distributions cannot be reliably predicted 
by assuming a constant factor between the {n — l)-jet and ra-jet cases. We describe the 
nontrivial structure found in the jet-production ratios for the vector-boson pt distributions 
using a simple model that incorporates leading-logarithmic behavior and suppression due 
to phase-space and parton-distribution-function effects. Related to the shape differences is 
a strong dependence of the jet-production ratios on the experimental cuts. In particular, 
there is about a 50% difference in the jet-production ratios for the CDF and DO setups, as 
shown in tables IVIII and IVIIII of section IVI[ 

This paper is organized as follows. In section [Til we briefly summarize our calculational 
setup, focusing on the differences from our previous work onW + 3-jet production [l7, 18|. 



Section UTTl records our choice of couplings, renormalization and factorization scales, and cuts 
matching those of the CDF and DO measurements. We also discuss issues associated with 
the choice of scale. In section IIV| we present results for cross sections and for a variety of 
distributions, matching CDF's cuts, and we compare with their measurements. In section IVl 
we give distributions in the softest jet pt for 7* + 1,2, 3-jet production, using DO's cuts 
and comparing with their published data. In section |Vl] we present jet production ratios for 
both the total cross section and the vector-boson pt distribution, and discuss a simple model 
for the latter ratio. We present our conclusions in section IVIII We include one appendix 
defining observables, as well as one giving values of the matrix elements at a single point 
in phase space. The latter appendix should aid future implementations of the Z, 7* + 3-jet 
one-loop matrix elements. 
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II. CALCULATIONAL SETUP 



A. Processes 

In this paper we calculate the inclusive processes, 

pp Z,7* + rajets + X ^ e+e" + n jets + X, (2.1) 

at ^/s = 1.96 TeV to NLO accuracy, for n = 1,2,3. Both CDF and DO have measured 

n n 

production cross sections for all three of these processes [2|, |^ at the Tevatron, based on 
integrated luminosities of 1.7 fb^^ and 1.0 fb"^, respectively. The DO measurements, besides 
using slightly different cuts, use a significantly smaller jet cone size, R = 0.5 versus R = 
0.7 for CDF. In addition, differential distributions have been provided for Z,'y* + 1,2- 
jet production. The set of available distributions is particularly extensive in the case of 
one jet ^, ^. The 2', 7* + 1, 2-iet production measurements have been compared to NLO 
predictions from MCFM 4l|. For the case of three additional jets, the data sets analyzed 
are still small, so that CDF has measured only a total cross section, while DO has provided 
three bins in the distribution of the transverse momentum of the third jet {pt ordered). 
The latter distribution was also compared with a leading-order prediction computed using 
MCFM. 

The present article provides the first NLO predictions for Z, 7* + 3-jet production, al- 
lowing a comparison with both the CDF cross section and the DO pt distribution. We will 
provide a few other distributions as well. We hope that, as additional data is analyzed, 
such distributions will be measured by CDF and DO. A comparison between NLO results 
and Tevatron data for various W, Z + 3-jet distributions would provide a very important 
benchmark for future LHC studies of these complex final states. 

In more detail, the process under consideration (12.11) receives contributions from several 
partonic subprocesses. At leading order, and in the virtual (one-loop) NLO contributions, 
these subprocesses are all obtained from 

qqQQg^Z,^* ^e+e- , (2.2) 
QQggg ^ ^, 7* ^ e+ e" , (2.3) 

by crossing three of the partons into the final state. The quarks are represented by q and Q 
and the gluons by g. The Z or photon couples to the quark line labeled q. Representative 
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FIG. 1: A few representative diagrams contributing to the qg — t- e^e~ qgg and qQ — t- e~^e~ qgQ 
one-loop amplitudes. The e~^e~ pair couples to the quarks via either a Z boson or an off-shell 
photon. 





g 




(b) 



FIG. 2: Sample diagrams illustrating one-loop contributions to Z, 7* -|- 3-jet production where the 
vector boson couples directly to a quark loop, via either (a) a vector coupling, or (b) an axial vector 
coupling. These contributions are quite small for the corresponding process with one parton less, 
and therefore are not included in our calculation. 

diagrams for the virtual contributions are shown in fig. [1] We include the decay of the 
vector boson {Z, 7*) into a lepton pair at the amplitude level. The photon is always off 
shell, and the Z boson can be as well. For the Z the lepton-pair invariant mass, Mee, follows 
a relativistic Breit-Wigner distribution whose width is determined by the Z decay rate Tz- 
We take the lepton decay products to be massless. Amplitudes containing identical quarks 
are generated by antisymmetrizing in the exchange of appropriate q and Q labels. 

The light quarks, u, d, c, s, b, are all treated as massless. We do not include contributions 
to the amplitudes from a real or virtual top quark. Nor do we include the pieces in which 
the vector boson couples directly to a quark loop through either a vector or axial coupling, 
as illustrated in fig. [2J In Z, 7* + 2-jet production these pieces affect the cross section by 
under 0.3%. We therefore expect the omission of these pieces to have a small effect on the 
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7* + 3-jet results presented here, well below the residual NLO uncertainties of 10-20%. 





Q 




-Qi 
Q2 



Q2 



FIG. 3: Representative real-emission diagrams for the eight-point tree-level amplitudes, qg 
e-^e~qggg, qg e~^e~qgQQ, and qq e+e" QxQxQiQi- 



Besides the loop amplitudes, we need tree amplitudes for real emission contributions. 
The relevant subprocesses are 



m9999 ^ ^, 7* ^ e+ , 
qqQQgg ^ Z, 7* e+ e" , 
qqQiQiQ2Q2 ^, 7* e+e" 



(2.4) 
(2.5) 
(2.6) 



where all the physical processes are obtained by crossing four of the partons into the final 
state. Representative tree diagrams for these contributions are given in fig. [3l 

To compute the NLO corrections we use BlackHat and SHERPA, essentially following 



the same calculational setup described in ref. [18!] for the W + 3-jet process. We therefore 
discuss our setup only briefly, pointing out the few differences with ref. 13] • 



B. Setup 

The virtual contributions are evaluated with BlackHat, which is based on the unitarity 
method 21|. One-loop amplitudes are expanded in terms of a basis set of scalar integrals 
composed of box, triangle and bubble integrals, plus a rational remainder. The coefficients 
of box integrals are obtained from quadruple cuts by solving the cut conditions [23]. Coeffi- 
cients of bubble and triangle functions are then obtained using a numerical implementation 
of Forde's approach In this implementation, BlackHat uses a procedure related to 
that of Ossola, Papadopoulos and Pittau [26^ to subtract box contributions when determin- 
ing triangle coefficients, and to subtract box and triangle contributions when determining 
bubble coefficients. The basis scalar integrals are evaluated numerically using their known 
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analytic expression s |65l| . To obtain rational terms, we have implemented both loop- level on- 
24j . |25| and a numerical version of the "massive continuation" approach due 



shell recursion 

to Badger 66|], which is related to the Z)- dimensional generalized unitarity 6?| approach of 



Giele, Kunszt and Melnikov 



29| . The numerical version involves subtracting the contribu- 



tions of higher-point cuts rather than taking large-mass limits. It is similar to the numerical 
version of Forde's method |27'] for four-dimensional unitarity cuts, which is described in 



ref. 



18|. In that paper we used on-shell recursion for the leading-color terms, where speed 



is at a premium. For the simplest helicity configurations, on-shell recursion is implemented 
analytically and the results stored for numerical evaluation. For subleading color terms the 
massive continuation method was used because it is presently more flexible. For production 
runs in the current study, we used the analytic formula obtained via on-shell recursion for 
the leading-color amplitudes, and the massive continuation method for the remaining terms. 

As discussed in ref. [l8|, for efficiency purposes it is useful to compute the leading- 
color parts of the virtual contributions separately from the numerically much smaller, but 
computationally more complicated, subleading-color contributions. We follow the same 
division of leading and subleading color as in ref. |l8|, except that here we assign the pieces 
proportional to the number of quark flavors (nj) to the leading-color contributions instead 
of the subleading-color ones. This has the effect of somewhat reducing the size of the 
(already very small) subleading-color contributions, helping to reduce the number of phase- 
space points at which they must be evaluated. We add the leading- and subleading-color 
contributions at the end of the calculation to obtain the complete color-summed result. 



We refer the reader to refs. 



22 



for detailed descriptions of the primitive amplitude 



decomposition that we used. Alternative organizations of color, within the context of the 



unitarity method, may be found in refs. 



34 



69|. 



An important issue is the numerical stability of the loop amplitudes. In fig. IU we illustrate 
the stability of the full-color virtual interference term (or squared matrix element), day, 
summed over colors and over all helicity configurations for the two subprocesses uu — )■ 
e^e^uug and uu — )• e^e^ggg. The horizontal axis of fig. H] shows the logarithmic error. 



(2.7) 



for each of the three components: 1/e and e°, where e = (4 — D)/2 is the dimensional 

regularization parameter. In this expression aW^ is the cross section computed by Black- 
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FIG. 4: The distribution of the relative error in the virtual cross section for the two subprocesses 
uu —7- e^e~uug and uu — )■ e^e~ ggg. The horizontal axis is the logarithm of the relative error ()2.7p 
between an evaluation by BlackHat, running in production mode, and a target expression evalu- 
ated using higher precision with at least 32 decimal digits (or up to 64 decimal digits for unstable 
points). The vertical axis shows the number of phase-space points out of 100,000 that have the 
corresponding error. The dashed (black) line shows the term; the solid (red) curve, the 1/e 
term; and the shaded (blue) curve, the finite (e'^) term. 



Hat as it normally operates for production runs (switching from 16 decimal digits to higher 
precision only when instabilities are detected), whereas Oy^^'^^ is a target value computed by 
BlackHat using multiprecision arithmetic with at least 32 digits, and 64 digits if the point 
is deemed unstable using the criteria described in refs. [l8|, l3l|- We use the QD package [70 1 



for higher-precision arithmetic. The phase-space points are selected in the same way as those 
used to compute cross sections. We note that an overwhelming majority (above 99.9%) of 
events are accurate to better than one part in 10^ — that is, to the left of the '—3' mark on 



the horizontal axis. Because we need only recompute parts of amplitudes in most cases 18|, 
the extra time spent in higher-precision operation is quite small, roughly 20% more than if 
only double precision had been used. 

In addition to the virtual corrections, the real-emission corrections are also required. 
These terms arise from tree-level amplitudes with one additional parton: an additional 
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gluon, or a quark-antiquark pair replacing a gluon, as illustrated in fig. [31 We use SHERPA 



for these pieces. The infrared singularities are canceled between rea^ 
contributions using the Catani-Seymour dipole subtraction method 



-emission and virtual 



42| , implemented 14 



in the automated program AMEGICH — h |^ 
We follow the same setup described in ref. 



|, which is part of the SHERPA framework |l3 |. 
18| . taking Odipoie = 0.03 as our default value. 



The Monte Carlo integration over phase space of both the real-emission and virtual pieces 



are carried out by SHERPA using a multichannel [7l| approach. In this approach, the in- 
tegrand is not split up into pieces, but is sampled differently in different channels. For 
Z, 7* + 1,2-jet production, we use AMEGIC-] — h, and each channel generates a momen- 
tum configuration based on the size of the denominators of the propagators of a tree-level 
Feynman diagram (Born or real-emission, as appropriate). For the more complicated case 
of Z, 7* + 3-jet production, in order to improve the efficiency, we have developed a spe- 
cific phase-space generator, applicable to + n-jet production. In this approach, a single 
channel generates a momentum configuration for the partons that is based on the size of de- 
nominators associated with a specific parton color ordering (the color-ordered QCD antenna 

nn 

radiation pattern), following the ideas of refs. [46|, |47[. The lepton momenta are generated 
so that the invariant mass of the lepton pair traces a Breit-Wigner distribution about the 
vector-boson mass. For larger numbers of partons this generator has a greatly reduced num- 
ber of channels, compared to the number of channels based on Feynman diagrams, so that 
it remains viable for vector-boson production with up to five or six jets. 

For the Z, 7* + 3-jet process we integrate the real-emission terms over about 4 x 10^ 
phase-space points, the leading-color virtual parts over 7 x 10^ phase-space points and 
the subleading-color virtual parts over 4 x 10^ phase-space points. The LO and dipole- 
subtraction terms are run separately with 10^ points each. These numbers are chosen to 
achieve an integration uncertainty of 0.5% or less in the total cross section. 

As a cross-check, we have compared our results for Z, 7* -|- 0, 1, 2-jet production at NLO 
and Z, 7*+ 3-jet production at LO to those of MCFM and find agreement to better than 1%. 
For Z, 7* + 2-jet production we used the same analytic one-loop matrix elements [22] as used 
in MCFM, with cross checks against a purely numerical computation within BlackHat. 
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parameter 


value 


aQED(Mz) 


1/128.802 


Mz 


91.1876 GeV 




0.230 




2.49 GeV 



TABLE I: Electroweak parameters used in this work. 
III. COUPLINGS, EXPERIMENTAL CUTS AND SCALE CHOICES 

In this section we describe the basic parameters used in this work, including couphngs, 
experimental cuts and our choice of renormalization and factorization scales. We also discuss 
the residual scale dependence remaining in the NLO results. 



A. Couplings and parton distributions 

We express the Z-boson couplings to fermions using the Standard Model input parameters 
shown in table HI The parameter g^, is derived from the others via, 



2 _ AnoQEDiMz) 

-2/1 

sm 6w 



9i= 27n ■ (3-1) 

We use the CTEQ6M 



72j parton distribution functions (PDFs) at NLO and the CTEQ6L1 
set at LO. The value of the strong coupling constant is fixed to agree with the CTEQ choices, 
so that as{Mz) = 0.118 and as{Mz) = 0.130 at NLO and LO respectively. We evolve as{fi) 
using the QCD beta function for five massless quark flavors for fi < rrit, and six flavors for 
fi > rrit. (The CTEQ6 PDFs use a five-flavor scheme for all fi > m^, but we use the SHERPA 
default of six-flavor running above the top-quark mass; the effect on the cross section is very 
small, on the order of 1% at larger scales.) At NLO we use two-loop running, and at LO, 
one-loop running. 

B. Experimental cuts for CDF 

To compare to CDF data we apply the same cuts as CDF j^, 
pJr>30GeV, |i/j"*|<2.1, 
14 



£;^>25GeV, |r/''|<l, (3.2) 

\r]^^\<l or 1.2 < |r/^2| < 2.8, 

Ai?e-jet > 0.7 , 66 GeV < M^e < 116 GeV . 

For any jet, pj^ denotes the transverse momentum and y^^^ the rapidity. For the leptons, 
denotes the transverse energy of either the electron or positron; 77^^ refers to the pseu- 
dorapidity of either the electron or positron and r/^^ refers to that of the other; M^e is the 
pair invariant mass. 



In their study of Z, 7* production, CDF used a midpoint jet algorithm 



6l| with a cone 



size of R = 0.7 and a merging /splitting fraction of / = 0.75. We use instead three different 



infrared-safe jet algorithms |55l-l57j : SISCONE (/ = 0.75), anti-Zcr and kx, all with R = 0.7. 
SISCONE is our default choice for comparison to CDF. (The fcy algorithm gives very similar 
parton- level results as the anti-Zcr algorithm, so we will not show those results explicitly.) 

Our calculation is a parton-level one, and does not include corrections due to non- 
perturbative effects, such as those induced by the underlying event, induced for example 
by multiple parton interactions, or by fragmentation and hadronization of the outgoing par- 
tons. In order to compare our parton-level results to data, we require non-perturbative 
correction factors. As discussed further in section [TV Bl for Z, j* + 1, 2-jet pt distributions, 
we adopt estimates of these correction factors made by CDF 



C. Experimental cuts for DO 



To compare to DO data we apply the jet cuts [4| 



pJ^*>20GeV, |r/J"'|<2.5 



,jet| 



(3.3) 



DO defined jets using the DO Run II midpoint jet algorithm 60|, with a cone size of i? = 0.5 
and a merging/splitting fraction of / = 0.5. We use instead the SISCoNE algorithm, with 
R = 0.5 and / = 0.5. 

DO performed an analysis with two distinct sets of lepton cuts. In their primary selection, 
which was compared directly to theory, only an invariant mass cut was imposed on the 
electron-positron pair. 



(a): 65 GeV < M^e < 115 GeV . 



(3.4) 
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For the secondary selection, the lepton cuts were, 

65 GeV < Mee < 115 GeV , 
(b): £;^>25GeV, (3.5) 
\r]^\ < 1.1 or 1.5 < \r]''\ < 2.5 . 

The latter "(b)" selection corresponds to the data DO actually collected. In their main 
selection ("(a)"), they extrapolated to an ideal detector with full lepton coverage using LO- 
matched parton-shower simulations. This extrapolation introduces an additional uncertainty 
and model dependence. It more than doubles the absolute cross section, although the 
quantities measured by DO, which are normalized by the inclusive Z, 7* + 0-jet cross section 
for the same lepton cuts, shift by much less. (Comparing the (a) entry to the corresponding 
(b) entry in table IIVI gives an estimate of the fraction of cross section in selection (a) that 
comes from the extrapolation.) 

We shall present NLO results corresponding to both selections, that is with and without 
the lepton acceptance cuts in the secondary selection (13.51) . Selection (b) allows us to 
compare to unextrapolated data. On the other hand, DO estimated the non-perturbative 
corrections, from hadronization and the underlying event, for selection (a) requiring 
us to extrapolate these corrections to selection (b) in order to use them there, as we shall 
discuss further in section |Vl 



D. Scale dependence 

Following the standard procedure, we test the stability of the perturbative results by 
varying the renormalization and factorization scales. In this article, we set the renormal- 
ization and factorization scales equal, = Hf = /W. In figs. [5] and E], we show the scale 
variation of the total cross section for the SISCONE and anti-fcr algorithms, respectively. 
In both cases we choose the central scale /xq = and then vary it down by a factor of 
four and upwards by a factor of eight. A fixed scale of order the Z mass is appropriate 
here, because the cross section is dominated by low-pj- jets. In both figures, the upper three 
panels show the markedly reduced scale dependence at NLO compared to the corresponding 
LO cross section in 7* + 1-, Z, 7* + 2-, and Z, 7* + 3-jet production, respectively. The 
bottom panel combines the ratios of NLO to LO predictions {K factors) for all three cases, 
illustrating the increasing sensitivity of the LO result with an increasing number of jets. 
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FIG. 5: The scale dependence of the LO (dashed blue) and NLO (solid black) cross sections for 
Z, 7* + 1,2,3-jet production at the Tevatron, as a function of the common renormalization and 
factorization scale with fiQ = Mz- Here the SISCone jet algorithm is used; the lepton and jet 
cuts match CDF [2]. The bottom panels show the K factor, or ratio between the NLO and LO 
result, for each of the three cases: 1 jet (dot-dashed red), 2 jets (dashed blue), and 3 jets (solid 
black) . 

This increase is expected, because there is an additional power of as for every additional jet. 
Accordingly, the reduction in the scale dependence at NLO tends to become more significant 
with an increasing number of jets. The plots for the kr algorithm are very similar to the 
ones for anti-/cT, so we do not show them here. 

Figs. |5] and |6] also reveal two further features. First, for n > 1 the cross section for the 
anti-/cT algorithm is significantly larger than for SISCoNE at the same value of R, especially 
at LO; the difference lessens at NLO. Second, the K factor at yU = Mz decreases significantly 
with the number of jets. The first feature is due to the smaller probability of two partons 
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FIG. 6: The scale dependence for Z, 7* + 1,2, 3-jet production at the Tevatron. The plot is the 
same as fig. [U except here the anti-/cT jet algorithm is used. 

clustering into a jet in the anti-fcr algorithm. In that algorithm (or in the algorithm), 
no clustering can take place unless the two partons are separated by less than R in the 
(?7, 0) plane; whereas in SISCONE they can be clustered out to a distance of 2R. Hence the 
effective area of a cone algorithm, for the same value of R, is somewhat larger (by a factor 



of about 1.35) than that of a cluster algorithm such as anti-A;T or kr 5J, |58|, |73|. At LO, 
clustering always causes a loss of events, and thus a decreased cross section for SIS Cone, 
relative to anti-^T- NLO corrections, however, tend to increase the cross section more for jet 
algorithms with larger effective cone areas, because there is less chance of radiating a parton 
out of the cone and thereby reducing the jet pt below the cut threshold 0, Q- Hence the 
cross-section difference between the algorithms is lessened at NLO. The differences between 
SIS Cone and krp algorithms at LO and NLO can also be examined as a function of the 



number of jets in + 1, 2, 3-jet production, using results for the LHC presented in ref. 



18| 



18 



However, in this work R = OA was used, resulting in far smaller perturbative differences 
between the algorithms. 

The second feature, that the K factors at /i = Mz decrease with the number of jets, 
is not unrelated. It was previously observed that W + 3-jet production for R = 0.4 had 
quite a small K factor 17|, [l8|, l36|]. The dependence on the number of jets was discussed in 
ref. [7^, where it was attributed to the LO cross section being "too high", in part because of 
collinear enhancements associated with the small jet size. We can see from figs. [5] and [6] that 
the trend of decreasing K factor is stronger for the anti-^T algorithm than for SISCONE. 
This feature is consistent with the picture of ref. Jt^, because the anti-Zcy effectively has a 
smaller jet size. 

For distributions, rather than total cross sections, we would like to choose a characteristic 
renormalization and factorization scale on an event-by-event basis, in particular to ensure 
that the tails of distributions are described properly. Previous studies (see e.g. refs. [2], 



16|, 



13, 



50( 1) have used the transverse energy of the vector b oso n, Erp as a common 



renormalization and factorization scale. As already argued in refs. 18|, ISlI] this choice is 
quite poor at LHC energies. Indeed, because of the large dynamic range at the LHC, at 
NLO the choice can go disastrously wrong for some distributions, leading to negative cross 
sections [1^. It also causes large changes in shapes of generic distributions between LO 
and NLO. This behavior reflects the emergence of large logarithms In fi/E, which spoil the 
validity of the perturbative expansion when n does not match the characteristic energy scale 
E. We note that without an NLO result for guidance, it may not be clear that a given scale 
choice — such as E^ — is problematic. 

Even for the Tevatron, with its smaller dynamic range, the commonly-used scale choice 
H = Ej^ is not particularly good. It leads to a large change in shape between LO and NLO 
in the px distribution of the third-hardest jet in Z, 7* + 3-jet production, as shown in the 
left panel of fig. [7l In contrast, the scale choice /i = Hrf^, where Ht is the total partonic 
transverse energy defined in eq. ( ]A4I) . results in little change in shape between LO and NLO. 
This choice is shown in the right panel of fig. [71 

The difficulty with using the vector-boson transverse energy, fi = E^, as the scale can be 
exposed jisj] by considering the two configurations depicted in fig. [HI In configuration (a), 
the Z boson has a transverse energy larger than that of the jets, and sets the scale for the 
process. In configuration (b), the two leading jets roughly balance in px, while the Z has 
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FIG. 7: The NLO px distribution of the third jet in 7* + 3-jet production at the Tevatron, for 
the SISCONE algorithm. In the left panel the scale choice // = is used and for the right panel 
/i = Ht/"^- The thin vertical lines (where visible) indicate the numerical integration uncertainties. 
The bottom part of each panel displays various ratios, where the denominator is always the NLO 
result at the reference scale choice, and the numerator is obtained by either evaluating the LO result 
at the same scale (dashed blue line), varying the LO scale by a factor of two in either direction 
(cross-hatched brown band), or varying the NLO scale in the same way (gray band). Although 
the two NLO results are compatible, the LO results have large shape differences, illustrating why 



/i = Ej: is not a good choice at LO. The jet and lepton cuts match those used by CDF 



much lower transverse energy. Here, the scale /i = is too low, and not characteristic 
of the process. In the tails of fig. [TJ configuration (b) dominates, because it results in 
a larger third-jet Pt for fixed center-of-mass partonic energy; contributions from higher 
center-of-mass energies which might boost the Z-boson transverse energy are suppressed by 
the fall-off of the parton distributions. This explains the large deviation between LO and 
NLO visible in the left panel of fig. [71 

In contrast, Ht (or some fixed fraction of it) does properly capture both configurations 
(a) and (b). It is thus a much better choice of scale. (For the purposes of fixing the scale, 
we prefer the partonic definition of the total transverse energy over the experimental one in 
eq. flA5l) because it is independent of the experimental cuts and the jet definitions [3].) In 
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Jl 

(a) (b) 

FIG. 8: Two distinct Z + 3 jet configurations with rather different values for the Z transverse 
energy. In configuration (a) an energetic Z balances the energy of the jets, while in (b) the Z is 
relatively soft. Configuration (b) generally dominates over (a) when the transverse energy of the 
third jet gets large. 

the remainder of this paper we take fi = Hrf^ as our default for both the renormalization 
and factorization scales, except where noted. To assess the remaining scale dependence in 
the cross sections we evaluate them at five scales: /i/2, /i/ a/2, /i, A/2/i, 2/i. We generate scale 



variation bands using the minimum and maximum values. In our previous analysis [18|] of W 
production in association with jets at the LHC, we chose /i = Ht- Generally, Ht tends to be 
on the high side of typical energy scales, so here we divide by a factor of two. The difference 
between the two choices at NLO is not large, on the order of 10% in the normalization, and 
with very small effects on the shapes of distributions. At LO the changes are, of course, 
larger, with up to 40% variations. 

Although we adopt here fi = Hrf^ as a good representative overall scale for general 
distributions, other approaches may be superior for particular distributions, or particular 
regions of phase space. For example, it may be possible to resum large logarithms that 
appear in particular corners of phase space, and match the resummed result to the NLO 
one. Even if that cannot be done, it is certainly possible that choosing a scale that is a blend 
of different scales (such as the different jet transverse momenta) is appropriate in some cases. 

In principle, perturbative approximations can break down in various kinematic regions, 
so it is important to check whether this can affect our results. The breakdown is often 
due to effects of soft-gluon emission that in can be resummed in many cases. Large soft- 
gluon effects can be obtained when there are explicit vetoes on soft radiation, or when such 
radiation is implicitly vetoed by fast-falling parton distribution functions. In this paper, we 
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put no explicit vetoes on soft radiation in any observable we consider. However, there is an 
implicit veto as one goes out in the tail of the Ht distribution or the third-jet pt distribution. 
This implicit veto might lead to large double Sudakov, or threshold, logarithms. 

In order to investigate whether such logarithms might be large, we take advantage of 
the fact that threshold logarithms in the high-p^- tail of the pt distribution for inclusive jet 
production should be very similar to the tails we are looking at in + n-jet production, at 
comparable values of jet pt- In both cases there is comparable mix of partonic channels, and 
similar values of parton x. Note, however, that one can reach much higher p^s experimentally 
in pure jet production because of the much larger cross sections. One recent resummation of 
threshold logarithms for inclusive jet production 52| shows that the effects are quite modest. 



For example, figure 6 of ref. 



52| shows the ratio K of the (matched) NLL result to the NLO 



result for single-inclusive jet production at the Tevatron Run I {^/s = 1.8 TeV), for pt from 
50 to 500 GeV. For various choices of the renormalization and factorization scales, K ranges 
from 0.98 to at most 1.14, as long as pt < 300 GeV. Note that 300 GeV is well above the 
third-jet »ts shown in fig. [71 (The relevant parton x values probed at the Tevatron in figure 



6 of ref. [5^ also correspond at the LHC to pts that are about seven times larger, well above 
the range studied in figure 9 of ref. [18]. There, the NLO cross section evaluated at /i = Ej^ 
became negative for a second-jet Et of only 475 GeV. Hence, even in this more extreme 
example, threshold logarithms are very unlikely to play a role in this behavior.) 

There is one other type of logarithm in ^ + n-jet production, which is not present in 
inclusive jet production, and that is a (double) logarithm of the form ln(pT',jct/^y), due to 
emission of electroweak vector bosons that are soft and coUinear with respect to the jets, as 
in the configuration shown in fig. [8](b). The importance of this logarithm was emphasized 
very recently 



53l | for the case of Z + 1-jet production. Although the NLO correction to this 
process is enhanced by In^^pTjet/My) with respect to LO, the effect is peculiar to + 1- 
jet production. It does not happen when two or more final-state partons are present at LO, 
because then the configuration in fig. Mh) can already be reached at LO. Also, because it is 
associated with electroweak boson emission, it does not represent a QCD double logarithm 
that will reappear at higher orders in a^. 

We conclude that there is no indication of a breakdown of fixed-order perturbation theory 



for the ranges of observables studied in this paper or in ref. 



18|. 
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IV. RESULTS FOR CDF 



In this section we present results for Z, 7* + 1, 2, 3-jet production (inclusive) at the Teva- 
tron, and compare to data from CDF that has been corrected back to hadron level j2|. For jet 
Pt distributions in Z, 7* + 1, 2-iet production, we use the (relatively large) non-perturbative 
corrections estimated by CDF 3] to transform our parton-level results to hadron-level ones. 
For Z, 7* + 3-jet production, non-perturbative corrections were not explicitly presented by 
CDF, and we give only parton-level results. 



A. Total cross sections 



# of jets 


CDF 
midpoint 


LO parton 
SISCONE 


NLO parton 
SISCONE 


LO parton 
anti-kT 


NLO parton 
anti-fcy 


1 


7003 ± 1461^^^ ± 406 


4635(2)1??« 


6080(12)1354 


4635(2)1??« 


5783(12)111 


2 


695 ± 37l^^ ± 40 


429.8(0.3)1JI|:I 


564(2)1^,9 


481.2(0.4)1}^^ 


567(2)131 


3 


60± lll^ ±3.5 


24.6(0.03)1^;^2^ 


36.8(0.2)1^1 


37.88(0.04)1^2-2 


44.7(0.24)1^;! 



TABLE II: 7* + 1,2, 3-jet production (inclusive) cross section (in fb) at CDF. The column 
labeled CDF gives the hadron-level results from ref. [2], using a midpoint jet algorithm. The 
experimental uncertainties are statistical, systematics (upper and lower) and luminosity. The 
columns labeled by LO parton and NLO parton contain the parton-level results for the SISCone 
and anti-ZcT jet algorithms. The central scale choice for the theoretical prediction is /i = Ht/2, the 
numerical integration uncertainty is in parentheses, and the scale dependence is quoted in super- 
and subscripts. Non-perturbative corrections should be accounted for prior to comparing the CDF 
measurement to parton-level NLO theory. 



In table HI] we present the total inclusive cross sections for 7* + 1,2, 3-jet production, 
showing both the CDF measurement and theoretical predictions, using our default scale 
choice /i = i^T/2. The theoretical results in the table are given for both the SISCoNE and 
anti-fc-r jet algorithms. (The kjn algorithm gives identical results as anti-A;^ at LO, and is 
within 1% at NLO.) In the second column we give the CDF measurement, for its midpoint 
jet algorithm and corrected to hadron level, along with the experimental uncertainties. 
The statistical, systematic (upper and lower) and luminosity uncertainties are given after 
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# of jets 


CDF 
midpoint 


LO parton 

SISCONE 


NLO parton 

SISCONE 


LO parton 
anti-fcy 


NLO parton 

anti-fcy 


1 


7003 ± Uett^l ± 406 


4206(2)1801 


6076(9)1^01 


4206(2)1801 


5828(9)1^25 


2 


695 ± 371:^1^ ± 40 


422.2(0.3)1}^^ 


576(2)1^? 


469.4(0.4)1}^^ 


583(2)1^} 


3 


60± llt^ ±3.5 


28.66(0.03)1}^;^ 


40.28(0.2)11^ 


43.3(0.05)126;6 


48.7(0.3)1^;^ 



TABLE IIL 7* + 1,2,3-jet production cross section (in fb) at CDF. This table is similar to 
table nil except that here the scale choice 



the central values. The third and fourth columns present the LO and NLO parton-level 
predictions. Here we quote the uncertainties from integration statistics in parentheses, and 
the scale dependence in super- and sub-scripts (upper and lower). The scale dependence 
is determined following the traditional prescription of varying the scale by a factor of two 
around the central choice /i = Ht/2, as described above. 

To assess the effect of changing the jet algorithm, we compare the SISCoNE and anti-kr 
results in table [Tll which further quantifies the differences that were visible in figs. [5] and 
[6l which used a fixed scale fi. Although the SIS Cone algorithm gives noticeably different 
results from the anti-fcy algorithm, the variations are similar in magnitude to the residual 
scale dependence. The reasons for the perturbative differences between SISCoNE and anti- 

algorithms were outlined in section IIII D[ 

It is also interesting to compare the results of table [TTl to cross sections obtained with the 
widely-used scale /i = instead of our default choice /i = Ht/'^- In table [Till we give cross 
sections with this scale choice for the SIS Cone and anti-fcr jet algorithms. Comparing 
these results to those of table Hit we see that, at least for 2', 7* +1,2 jets, the K factor 
(ratio of NLO to LO) is much closer to unity for the choice /i = than for /i = E^. 

Although the /i = E^ choice is problematic in general, as already noted in section IIIIl at 
NLO it gives results for the total cross section that are similar to those from our default 
choice of /i = Ht/2. 

In order to compare parton-level results to the experimental measurement we must ac- 
count for non-perturbative corrections, using estimates from CDF [2] . These corrections are 
sizable, increasing the total cross section by a factor between 1.1 and 1.4 as the number of 
jets increases from one to three. As we will see below in fig. [9] for the jet pt distributions, 
these estimated correction factors align NLO theory with the measurement within uncer- 
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tainties, although a much more careful study of the non-perturbative corrections and the 
differences in jet algorithms is needed. It is interesting to note that in the CDF measurement 
of + n-jet production [l^, the corrections are significantl y sm aller. That measurement 
used a jet cone size of i? = 0.4 (with the JETCLU algorithm [75]). There the hadronization 
and underlying event corrections were under 10% below 50 GeV and under 5% at higher 
Et- The CDF study may also be contrasted with the DO study {4] discussed below, in 
which the cone size of R = 0.5 leads to non-perturbative corrections on the order of 15%. 
From the perspective of maintaining the precision of NLO predictions, it is advantageous 
to choose jet-cone sizes which minimize non-perturbative corrections, while not increasing 
the size of (In /^-enhanced) higher-order perturbative corrections too much. As discussed in 



e.g. refs. 
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73| . there is a tradeoff between the underlying event correction (increases as 



R increases) and splash-out (increases as R decreases), and a careful study would be needed 
to find the best choice. 



B. Comparison to CDF jet px distributions 

In this subsection, we compare our results with CDF data for jet pt distributions in 
Z, 7* -|- 1-jet and Z, 7* + 2-jet production. In the observables used by CDF, sometimes 
referred to as inclusive jet px distributions, all jets passing the cuts are included in the 
distributions. That is, if n jets pass the cuts, the event is counted n times, with contributions 
to each of the n bins containing the transverse energy of one of the jets. By definition, for 
inclusive Z + n-jet production, at least n jets pass the cuts, and periodically additional jets 
can also pass the cuts. At NLO these extra jets are modeled by a single extra jet from the 
real-emission contribution. This causes the area under the curve to be slightly more than n 
times the total cross section. In contrast, the W + n-jet production distributions measured 
in ref. [l^ and the Z + n-jet production distributions measured in ref. J] are differential in 
the transverse energy (or momentum) of the n*^ jet, and each event is counted only once, 
so they integrate to give the total cross section for V + n-jet production. 

The left and right panels of fig. [9] are for Z, 7* + 1-jet and Z, 7* + 2-jet production, 
respectively. The upper part of each panel compares the LO and NLO results against CDF 
data from ref. ^]. BlackHat+SHERPA produces NLO parton-level predictions. To 
compare to the CDF measurement we need to account for non-perturbative corrections. We 
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FIG. 9: Jet pT distributions for 7* + 1,2-jet production at the Tevatron with CDF's cuts. 
The theoretical predictions use the SISCone algorithm and the scale choice /i = Ht/2. In the 
upper panels the parton-level NLO distribution are the solid (black) histograms, and the NLO 
distributions corrected to hadron level are given by dash-dot (magenta) curves. The CDF data 
are the (red) points, whose inner and outer error bars denote respectively the statistical and 
total uncertainties on the measurements (the latter obtained by adding separate uncertainties in 
quadrature). The LO predictions have been corrected to hadron level and are shown as dashed 
(blue) lines. The lower panels show the distributions normalized to the full hadron-level NLO 
prediction for /x = Ht/2. The scale-dependence bands in the lower panels are shaded (gray) for 
the NLO prediction corrected to hadron level and cross-hatched (brown) for LO corrected to hadron 
level. 

use the last column of Table I of ref. l2|l as an estimate of their size. This table of corrections 

n 

was determined for the CDF midpoint jet algorithm using Pythia an LO-based parton- 
shower, hadronization and underlying-event program. Because we used the (infrared-safe) 
SISCone algorithm, the possible algorithm-dependence of the non-perturbative corrections 
introduces additional uncertainty into the comparison. As mentioned in the introduction, 



studies 
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62| of inclusive-jet cross section differences between midpoint algorithms and 
SISCone (which were also performed for R = 0.7) suggest relatively small "parton-level" 
differences between the algorithms, which in turn suggests that applying the CDF non- 
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perturbative corrections to our SISCoNE perturbative prediction is not unreasonable. The 
size of the corrections can be seen in the upper panels of fig. [9] by comparing the curves 
labeled "NLO parton", which are the parton-level predictions, to the ones labeled "NLO 
hadron", which are the hadron-level ones. It is easier to judge the size in the lower panels, 
using the solid (black) curves which give the ratios of the two predictions. For example, 
for Z, 7* + 2-jet production, non-perturbative corrections are significant for low px-, on the 
order of 20% at 30 GeV, and gradually drop to under 5% at larger jet transverse momenta. 
Uncertainties in the non-perturbative corrections are not included in the plots. 

The bottom panels shows various ratios, normalized to the NLO hadron-level prediction 
for the central scale /i = -ffr/2. We include scale-dependence bands, as described above, 
for the predictions corrected to hadron level. As expected, for NLO the scale dependence is 
greatly reduced when compared to LO. We note that for both 7* + 1, 2-jet production, 
the NLO hadron-level jet pt distributions match the CDF results quite well, noticeably 
better than the hadron-level LO distributions or parton-level NLO distributions. A similar 
comparison of the experimental data to NLO predictions was given in the CDF study, using 
MCFM [41']. The ratios of data to NLO presented there differ by up to 10% from those 
shown in fig. [9l Most of the difference can be attributed to the choice of central scale in the 
NLO result, ^ = versus /i = Ht/'^- CDF also assessed the uncertainties on the NLO 
predictions arising from the parton distribution functions. They found them to vary from 
4% at low jet pt to 10% at high pr, which is generally smaller than the NLO scale variation. 

C. Predictions for Z, 7* + 3-jet distributions at CDF 

In fig. [TOl we show the combined distribution of all jet p^s in Z, 7* + 3-jet production. 
It would be very interesting to compare this prediction to CDF data, after accounting for 
non-perturbative effects. As discussed above, the integral under the curve gives a bit more 
than three times the total cross section. As can be seen in the plot, with the scale choice 
= Ht/'^ there is only a modest change in shape between LO and NLO, especially at higher 
jet Pt- This is similar to the parton-level results at Z, 7* -|- 1, 2-jet production shown in fig. [91 
We expect non-perturbative corrections to lead to larger shape changes at lower pt- 

The separate distributions for the hardest, second-hardest, and third-hardest jet are 
shown in fig. [TTJ The shapes of the LO and NLO distributions are again similar, with our 
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FIG. 10: pt of all jets for 7* + 3-jet production with the CDF setup, using the SISCone jet 
algorithm and the scale choice ^ = Ht/2, for LO and NLO at parton level. The thin vertical 
lines (where visible) indicate the numerical integration uncertainties. The lower panel bands are 
normalized to the central NLO prediction, as in fig. [9j 



default scale choice. As in W+ 3-jet production [18|, successive jets have increasingly steeply 
falling distributions. 

In fig. [12] we show the rj distribution of the positron for .Z", 7* + 1,2, 3-jet production. 
The discontinuity and gap between 77 = ±1 and r/ = ±1.2 result from the discontinuity in 
the charged-lepton cuts in eq. (13. 2p . A careful inspection reveals a small forward-backward 
asymmetry, which can be traced to the left-right asymmetry in the Z boson couplings to 
fermions. (Similar asymmetries are discussed in e.g. ref. Once again, there is only a 

modest shape change between LO and NLO. 

V. RESULTS FOR DO 



The DO collaboration has studied jet px distributions in inclusive Z, 7* + n-jet production 
for up to three jets [^], using the DO Run II midpoint jet algorithm. Here we present 
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FIG. 11: First, second and third jet pT distributions for Z, 7* + 3-jet production. The dashed 
(blue) hnes are LO predictions and the solid (black) lines are NLO predictions. The SISCone jet 
algorithm and a scale choice of /x = Ht/'^ are used for these plots. 



# of jets 


LO parton 
selection (a) 


NLO parton 
selection (a) 


LO parton 
selection (b) 


NLO parton 
selection (b) 


[pb] 


179.01(0.02)1°:649 


236.96(0.08)1|J1 


84.08(0.03)1°:^! 


106.81(0.15)1°:!^ 


1 [fb] 


25223(9)11011 


30230(55)1^^^? 


10083(6)1}927 


12537(44)1^^° 


2 [fb] 


3787(3)11^9 


4415(14)1^60 


1538(2)16°8 


1848(10)1127 


3[fb] 


462(0)11° 


553(3)1™ 


190(0)lif 


236(1)1!° 



TABLE IV: NLO parton-level Z,^* + 0,1, 2, 3-jet production cross sections corresponding to DO 
selections (a) and (b). The columns labeled by LO parton and NLO parton correspond to the 
parton-level results for the SISCone algorithm. The central scale used for one or more jets is 
= Ht/"^- The numerical integration uncertainties are in parentheses and the scale dependence 
in super- and subscripts. 

the corresponding NLO parton-level results. To compare NLO theory and experiment we 
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FIG. 12: The r] distribution of the positron, for Z, 7* + 1,2,3-jet production. The discontinuities 
in the plots are due to the experimental cuts (j3.2p . The cross sections in the second and third 
panels are multiplied by factors of 10 and 100, respectively. 

again need to account for non-perturbative corrections. DO has provided estimates of non- 
perturbative corrections due to the underlying event and hadronization effects for their 
study using the lepton cuts (a) of eq. fl3.4l) . With the smaller cone size used by DO, R = 0.5, 
the net non-perturbative corrections turn out to be no larger than about 15%, significantly 
smaller than in the CDF analysis for R = 0.7. As described further below, we will use these 
correction factors as a rough estimate of the non-perturbative corrections for selection (b) 
as well. 
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A. Total cross section 



The theoretical predictions for the total cross sections for selection (a), with lepton 
cuts (13 ■4p . and for selection (b), with lepton cuts fl3.5l) . are given in table [IVl The LO 
and NLO parton-level cross-sections are for the SISCONE algorithm, with the central scale 
choice /i = and the scale dependence determined as before. For the case of Z, 7* + 0- 

jet production we use fi = E^, because /i = Hrf^ can vanish. As seen from table HVl for 
Z,j* + 1,2,3-jet production the LO scale dependence is quite large, but is substantially 
reduced at NLO. In particular, for Z, 7* + 3-jet production a shift in the scale by a factor of 

two causes a variation of up to 60% at LO and under 20% at NLO. For the case of Z, 7* + 0- 

f — i 

jet production, the scale dependence of the LO result is anomalously small [77|. This scale 
independence is due to the absence of factors of the strong coupling as in the LO matrix 
elements, along with the mild dependence of the quark distribution functions q^Xjfip) on 
/ii? at values of x that are relevant at the Tevatron.^ 

In ref. DO provided jet px distributions rather than summed cross sections. We could, 
of course, use the differential measurement to obtain the total cross section for Z, 7* + 1,2, 3- 
jet production (as a fraction of the inclusive Z-boson production cross section), but as the 
systematic error correlations were not specified, we instead turn to a direct comparison of 
theory and experiment for the pt distributions. 



B. Comparison to DO jet px distributions 

We now compare our predictions for jet pt distributions to the DO measurements. From 
our vantage point, selection (a) allows a direct comparison of our NLO predictions to the 



plots presented in ref. {4], which also include results from various LO and parton-shower 
based programs. On the other hand, as noted earlier, in selection (a) more than half the 
events are extrapolated from the actual measurement using the cuts (13.51) of selection (b). 
This extrapolation introduces additional uncertainty into the measurement. Thus selec- 
tion (b) would normally be preferred for comparison to theory. However, because the non- 



^ The smaller scale variation of inclusive Z, 7* production at LO, with respect to NLO, does not, of course, 
imply that the LO prediction is more accurate; indeed, the NLO result is much closer to the NNLO 



one 



77|- 
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perturbative corrections determined in ref. j4| were for selection (a), they do not directly 
apply to selection (b). Selection (b) corresponds to a subset of (a), which could have some- 
what different average values for the Z boson rapidity and transverse momenta; these values 
could in turn affect the jet kinematics and therefore the non-perturbative corrections. Be- 
cause the correlation between lepton cuts and jet kinematics is at second order, and because 
the quoted non-perturbative corrections for (a) are no larger than about 15%, one may hope 
that the corrections for (b) are not too different. For either selection, we face the same 
issue as with CDF, that the non-perturbative corrections were estimated for a midpoint 
jet algorithm rather than SISCONE. In light of these issues, we present NLO parton-level 
results with the SISCONE algorithm as our primary predictions. We then adopt the non- 
perturbative corrections given in ref. [^J for both selections (a) and (b), leaving possible 
(significant) improvements to future studies. 

In the comparison, we follow DO in normalizing the 2", 7* + 1,2,3-jet distributions 
by dividing by the inclusive Z, 7* + 0-jet cross section. The latter cross section is defined 
using the same set of lepton cuts as applied in the numerator, for both the (a) and (b) 
selections. When assessing the scale dependence, we vary the scale by a factor of two in 
each direction in the numerator, but in the denominator for simplicity we always take the 
central-scale values in the first row of table IIV[ Because both the LO and NLO inclusive 
Z, 7* cross sections vary by under 2% from their central values, this procedure modifies the 
scale variation band only slightly, with respect to varying scales in the denominator as well. 

We have compared DO data for Z', 7* -|- 1,2,3-jet production to our NLO prediction for 
both selections (a) and (b). Ref. Q] already showed a comparison between selection (a) for 
inclusive 2', 7* -|- 1,2-jet production and NLO predictions using MCFM 4l|], so for these 
processes, we show instead a comparison to selection (b). The comparison is displayed in 
fig. [131 As explained above, we use tables IV and V of ref. to convert our parton-level 
results to hadron- level ones. Ref. ^| studied the effects of parton distribution function 
uncertainties on the perturbative predictions for selection (a), and found them to be 5-10% 
for the leading two jets, and 5-15% for the third jet. We have compared our results for both 
selections (a) and (b) to MCFM, using the fc-r algorithm with scale choice /i = we find 
agreement, with the total cross section agreeing to better than 0.5%. (As was the case for 
the CDF study, the ratio of data to NLO presented by DO differs by up to 10% from our 
corresponding results, though for selection (b), in fig. [T31 again the difference is largely due 
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FIG. 13: Normalized jet pT distributions for DO selection (b). The left plot shows l/az-f* x da/dpT 
for the leading jet in Z, 7* + 1-jet production. The right plot shows the distribution for the second- 
hardest jet in Z^^* + 2-jet production. The SISCone algorithm and scale choice /i = Ht/"^ are 
used in the theoretical predictions. In the upper panels the parton-level NLO distributions are the 
solid (black) histograms, while the NLO distributions corrected to hadron level are given by dash- 
dot (magenta) histograms. The DO data is indicated by the (red) points; the inner and outer error 
bars denote respectively the statistical and total experimental uncertainties on the measurements. 
The LO predictions corrected to hadron level are shown as dashed (blue) lines. Each lower panel 
shows the distribution normalized to the full hadron-level NLO prediction. The scale-dependence 
bands in the lower panels are shaded (gray) for NLO and cross-hatched (brown) for LO. 

to the different choice of central scale in the NLO result, = versus = Ht/2.) 

Fig. HH compares our theoretical predictions for the third-jet pt distribution in Z, 7* + 3- 
jet production to both selections (a) and (b). As expected, the scale dependence of the 
NLO predictions is quite a bit smaller than for LO, throughout the distribution, with only 
a 15% deviation from the central value. For reference, we also present the results shown in 
fig- EH in tabular form, in tables |V] and I VII The columns labeled "LO parton" and "NLO 
parton" give the parton-level results. These are the primary results of our DO study, and 
would be the key input to any future analyses. The columns labeled "LO/NLO hadron" 
give these predictions multiplied by the non-perturbative corrections given in table VI of 
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FIG. 14: Comparison of NLO theory to the DO resuh for the distribution of l/az,'y* x da/dpx for 
the third-hardest jet in inclusive Z, 7* + 3-jet production. The left plot shows the comparison for 
selection (a), and the right plot for selection (b). The quoted 3] non-perturbative corrections for 
selection (a) were used as estimated for both selections. The SISCone jet algorithm and the scale 
= Ht/"^ were used in the theoretical predictions. The labeling is as in fig. [T3j 
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TABLE V: Z,^* + 3-jet production for DO selection (a), l/az,^- ^ da/dpx for the third jet {pT 
ordered) in inclusive Z,j* ± 3-jet production, with the lepton cuts of eq. p.4p . The LO and 
NLO hadron columns include hadronization and underlying-event corrections from ref. B. The 
hadron-level columns do not include the uncertainties from non-perturbative corrections; estimates 



of these may be found in table VI of ref. 



a 



ref. and represent a rough estimate. As in fig. [Tl]we have not included the uncertainties 
in the non-perturbative corrections. Finally, the column labeled "DO" gives the DO measure- 
ment, followed by its statistical and systematic uncertainties. Fig. [T4] and tables IVl and NH 
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TABLE VI: 7* + 3-jet production, l/crz^-y* x da/dpT for the third jet {pT ordered) in inclusive 
7* ± 3-jet production. The setup is the same as in table |Vl except the additional lepton cuts 
given in eq. ()3.5p for selection (b) are imposed. In the columns labeled by "hadron" we again 
multiplied by the non-perturbative corrections (computed for (a)) from Table VI of ref. J]. 

show that, although the experimental central values are always a bit above the theoretical 
bands, the agreement between theory and experiment is reasonable, given the experimental 
statistical uncertainties, and despite the various unquantified uncertainties discussed above. 
For both theory and experiment, the shift in values between selections (a) and (b) is quite 
small, significantly smaller than the respective uncertainties. 

VI. JET-PRODUCTION RATIOS 



The measurement of jet cross sections is sensitive to a number of experimental and non- 
perturbative issues, in particular the measurement of the jet energy and contributions due 
to the underlying event. The latter is not modeled in perturbative predictions, and accord- 
ingly introduces a systematic uncertainty. Jet pr distributions fall rapidly: the distribution 
of the third-hardest jet in Z + 3-jet production in fig. [11], for example, falls by two orders 
of magnitude over a factor-of-two range in px- Thus small errors in pt measurements or 
non-perturbative shifts in jet Pt can have important effects on distributions. Furthermore, 
the cross section for an additional jet is roughly an order of magnitude smaller, so a misiden- 
tification of an {n — l)-jet process as an n-jet process can cause a significant error in the 
measured cross section for the process with the extra jet. These difficulties increase as the 
number of jets accompanying a vector boson grows. 



A simple way to contro 



of cross sections 
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some of the systematic uncertainties is to consider instead ratios 



78| . As an example, consider the jet-production ratio (also known as 



the "Berends" or "staircase" ratio), the ratio of Z, 7* ± n- to Z, 7* ± (n — l)-jet production. 
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TABLE VII: The ratios of the Z,^* + n-jet to Z,j* + (n — l)-jet cross sections for CDF's cuts, 
using the SISCone algorithm and /x = Ht/2. For the experimental uncertainties we removed the 
luminosity errors, as these cancel; we treated all remaining uncertainties as uncorrelated. For the 
theoretical ratios we varied the scale in the same way in the numerator and the denominator. The 
integration uncertainties are small compared to the remaining scale dependence. 
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TABLE VIII: The ratios of the Z,'j* ± n-jet to Z,j* ± (n — l)-jet cross sections for DO's selection 
(b), using the SISCone algorithm and fi = Ht/2. Here we keep only the experimental statistical 
uncertainties and drop the systematic ones. The scale dependence is treated as in the previous 
table. 

displayed in table IVTl] for the CDF cuts and in table IVml for DO type (b) cuts fl33|l . 
We expect such ratios to be much less sensitive to experimental systematic measurement 
uncertainties than the individual cross sections.^ Ratios also mitigate various theoretical 
uncertainties, such as uncertainties in the non-perturbative corrections. The theoretical 
values in table IVIII and table I VI iTl are given at parton level. In assessing the scale dependence 
of the ratio, we varied the renormalization and factorization scales (by a factor of two) in the 
same way for both numerator and denominator. (Of course, one could consider alternative 
schemes to estimate scale dependence.) Lacking knowledge of the correlations between the 
experimental systematic uncertainties, we took them to be uncorrelated. (We also took the 
larger of the upper and lower uncertainties, as they are quite close for all cases.) Both 

^ Note that the jet-energy scale uncertainty will not completely cancel, however, because the px distribution 
in the n"^ jet is steeper than that in the {n — 1)^* jet. 
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the NLO and LO results are compatible with the experimental ratios, within the estimated 
uncertainties. It would be interesting to re-examine these ratios using the larger data set 
collected more recently at the Tevatron, and properly incorporating all correlations in the 
experimental systematic uncertainties. 

Previous studies have noted that the jet-production ratios are roughly independent of the 
base number of jets. With the choice of cuts used by CDF, eq. (13.21) . our NLO theoretical 
results for n = 2,3 are similar, within 30% of each other. With DO's choice of cuts, the 
agreement is better, to within 15%. The predicted ratio for n = 4 is significantly lower, 
but this prediction is still only at LO, because the NLO results for 2', 7* + 4-jet production 
are not yet available. Tables IVIII and IVIIII also show that the NLO predictions for the 
jet-production ratios are quite close to the LO central values, but tend to have less scale 
variation. 
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FIG. 15: The LO and NLO vector-boson pT distributions for 7* + n-jet production at CDF. In 
the upper panels, the top distribution is for one-jet inclusive production, the one underneath it is 
for two-jet inclusive production, and the next one is for three-jet inclusive production. At LO, in 
the left plot the bottom curve is for four-jet production. The lower panel gives the jet-production 
ratios as a function of vector-boson pj-- In the bottom panel the upper curve is the 2/1-jet ratio 
and the one underneath it the 3/2-jet ratio. The 4/3-jet ratio at LO is displayed as well (magenta 
dotted curve). 
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FIG. 16: Comparison of jet-production ratios, differential in the vector-boson px, at LO and NLO. 
The lower panel shows these ratios, divided by the NLO ratio evaluated at the default central scale 
choice, and including scale variation bands (computed as for the total cross section ratios). 

This approximate independence, however, hides a great deal of variation in differential 
distributions. In fig. [151 we show the differential distributions in the vector-boson transverse 
momentum (p^) for inclusive 2', 7* + 1,2,3,4-jet production at LO (left panel) and for 
inclusive 2', 7* + 1,2,3-jet production at NLO (right panel). The lower panes show the 
corresponding ratios: (Z, 7*+ 2-jet)/(Z, 7*+ 1-jet) (2/1), (Z, 7*+ 3-jet)/(Z, 7*+ 2-jet) (3/2), 
{Z,Y + 4-jet)/(Z,7* + 3-jet) (4/3) at LO; 2/1 and 3/2 at NLO. These ratios have shapes 
that are quite stable in going from LO to NLO, as shown in fig. [161 with the exception of the 
2/1-jet ratio at low pt and small differences in the 3/2 ratio at higher pt- In these plots, we 
use parton-level results, without any corrections for hadronization or the underlying event. 
We expect substantial cancellations of these non-perturbative corrections in the ratios. 

This figure shows that the jet-production ratios depend strongly on the px of the vector 
boson, and that the 3/2-jet and 2/1-jet ratios are rather different beyond low pt- This 
means that their putative independence of the base number of jets is illusory: in reality, 
they depend sensitively on the cuts applied. For example, a p-r > 70 GeV cut on the vector- 
boson transverse momentum would result in a rather sizable difference between the 3/2-jet 
and 2/1-jet ratios of total cross sections. 
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How does the non-trivial dependence of these ratios on pj, arise? At LO, the 2/1-jet ratio 
is undefined (infinite) at low ; it rises smoothly from a very small value at — p^miin 
where Pr'min minimum jet set by the experimental cut (30 GeV in this case). In 

contrast, at NLO the 2/1-jct ratio rises to a finite value as — ?> 0. The other ratios have no 
structure at low px- All of the ratios rise noticeably beyond p}^ of 70 GcV or so, a rise which 
continues up to around 200 GeV, where the ratios fiatten out or start to decline somewhat. 

The behavior of the ratios at low pt is relatively easy to understand. For p"^ below the 
minimal jet pt, the vector's transverse momentum cannot be balanced by a lone parton, so 
the leading-order contribution to 7* -|- 1-jet production vanishes. Accordingly, the 2/1 
ratio is infinite at LO for pj. < Pxmin) ^^^^^ rises from a very small value just above p}f = PTmm- 
For Z, 7* + 1-jet production at NLO, as pj' 0, the only contribution is from real-emission 
configurations with two partons which are roughly balanced in transverse momentum. It 
is reasonably hkely that the second hardest parton also has p^ > p^^Lin- Therefore the 
differential cross section in this region is of the same order in as, 0{al), as the leading 
contribution to 7* -|- 2-jet production. Hence the NLO 2/1 ratio rises as — > 0, to a 
number independent of as (and "of order unity"). No such kinematic constraints arise in 
vector-boson production accompanied by more than one jet, even at LO, so the 3/2 and 4/3 
ratios remain small as p^ ^ 0. 

What about the rise at larger vector-boson transverse momentum? For a given large p^, 
we expect that the matrix element is maximized for an asymmetric configuration of jets, 
corresponding to a near-singular configuration of the partons. A typical configuration, for 
example, would have one hard jet recoiling against the vector, and additional jets (if any) 
with transverse momenta down near the pt cut. In these configurations, the short-distance 
matrix element will factorize into a matrix element for production of one hard gluon, and a 
singular factor: either a splitting function in coUinear limits, or an eikonal one in soft limits. 
The phase-space integrals over these near-singular configurations give rise to potentially 
large logarithms. Because the minimum AR is relatively large, coUinear logarithms should 
not play an important role; on the other hand, Pr/Prmm can become large, so its logarithm 
will play a role. 

The approximate factorization suggests that we can model the differential cross sections 
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shown in fig. [15] by the following forms, 



(74 



a^(6o + &ilnp)/(p^)(l-p^/pr'') 



\12 



aiico + ci Inp + C2 In^ p)/(p^)(l - p^/p^ 



V / max\73 



(6.1) 



atido + d,\np + d2 \n' p + ds \n' p)f{p^){l - p^/p"? 



V / max\74 



where a, = as{p^) N J {27c), p = {pX^/p^Tmin? and p??^^ = 980 GeV. The additional factors 
of (1 — p'j, / p^'^^)"' take into account the limits of the different-dimension phase spaces and 
possibly different parton-distribution-function suppression in the four cases. The function 
fiPr), which describes the overall rapidly-falling form of the distribution, will cancel in the 
ratios, leaving us with three parameters for the 2/1 ratio, four additional ones for the 3/2 
ratio, and five further parameters for the 4/3 ratio. 
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FIG. 17: Fits of the LO jet-production ratios (using CDF cuts) as a function of pT to forms derived 
from eq. ()6.ip . Top left, the 2/1 ratio; top right, the 3/2 ratio; bottom, the 4/3 ratio. The vertical 
scales are different. 
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We can determine these parameters for the CDF cuts by fitting the ratios to our LO 
results (where much smaller statistical uncertainties are easier to achieve). We evaluate the 
LO distributions at the central fj, value, and fit in the range 100 GeV < — 500 GeV. The 
result is 

60 = -1.642, 61 = 2.437, 72 = 1.118, 

Co = 5.081, ci = -5.812, cs = 2.658 , 73 = 2.539, (6.2) 
c/o = -5.728, c/i = 10.945, d2 = -6A, 4 = 1.628 , 74 = 4.195. 

The separate fits to the LO ratios are shown in fig. [T71 In spite of the limited number of 
parameters, the model gives an excellent approximation to the LO ratios, even down to 
of 70 GeV for the 2/1 ratio. (For the 2/1 ratio, the model's predictions are within 0.7% of 
the numerical results for two-thirds of the points in the fit range, and all points are within 
twice the Monte Carlo statistical error. For the 3/2 ratio, two-thirds of the points in the fit 
range are within 1%.) Although we did not attempt a fit to the NLO results for the 2/1 
or 3/2 ratios, due to larger statistical uncertainties, fig. [16] indicates that the parameters 
entering the 2/1 ratio would be essentially unchanged, and the 3/2 ratio parameters would 
only change modestly. We stress that the model is a purely phenomenological one, with 
some physically motivated input as to its form. It is not intended to supply a perfect fit, 
but it is remarkable that it does so well with so few parameters. 

The overall conclusion of this section is that the jet-production ratio is not solely a 
measure of as, but also depends sensitively on kinematical ratios. 



VII. CONCLUSIONS 



In this paper we have presented the first full NLO results for Z, 7* + 3-jet production at 
a hadron collider. (Ref. 48| contains a preliminary version of some of these results, based 



on a leading-color approximation.) We chose to present results for the Tevatron, so that we 
could compare to existing data from CDF and DO 2, J]. Our study should also serve as a 
benchmark for future NLO studies of Z boson production in association with multiple jets 
at the LHC. 

In line with expectations, we found that NLO results for Z, 7* + 1, 2, 3-jet production are 
much less sensitive than their LO counterparts to the choice of a common renormalization 
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and factorization scale /x. The deviation in the cross section, when varying around our 
default central scale choice by a factor of two, drops from 60% at LO to 15-22% at NLO. 



As noted in earlier studies 
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5l| . the commonly-used scale choice of the vector-boson 



transverse energy is not particularly good at LHC energies, as it leads to large shape changes 
in distributions between LO and NLO. At LHC energies, this choice can even lead to negative 
cross sections at NLO, in tails of various distributions [18]. As we have shown in the present 
paper, this choice is also a poor one for higher-multiplicity Z, 7* + n-jet production, even 
at Tevatron energies, because it generically results in large shape changes from LO to NLO. 
We used half the total transverse energy, /i = Ht/'^i as our default scale choice instead. 

The NLO results presented in this paper are parton-level results. To compare to the 
measurements, non-perturbative corrections need to be taken into account. We used the 
non-perturbative corrections tabulated by CDF and DO to estimate the corrections to our 
parton-level predictions. These corrections were computed for different cone algorithms, 
and in the case of DO selection (b), for a slightly different setup. We appealed to inclusive- 
jet studies [55I, |62| of algorithm differences to suggest that applying these corrections is 
reasonable; clearly, further study of this issue is desirable. Applying the corrections, we 
found good agreement between the NLO results and the CDF and DO measurements of 
cross sections and jet px distributions. For Z, 7* + 3-jet production, both CDF and DO data 
lie somewhat higher than the theoretical scale-dependence band, although the statistical 
uncertainties are still large. Obviously, from the perspective of comparing to NLO results, 
it is helpful to choose infrared-safe jet algorithms, as well as cone sizes for which the non- 
perturbative corrections and their associated uncertainties are small. 

For the cuts used by CDF, we presented results for two algorithms, SISCONE and anti- 
kx- The differences in cross sections for the two alg orithms can be explained qualitatively 



by the larger effective cone size of SISCoNE 
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731], for a given cone-size parameter R. 



We presented new results on jet-production ratios. We confirmed that for total cross 
sections the ratio of n-jet production to {n — l)-jet production is roughly constant 63|, l64 ]. 
For the cuts used by CDF, the NLO prediction gives ratios for n = 2, 3 that are within about 
30% of each other, whereas for the cuts used by DO the corresponding ratios are within about 
15%. (We note that the ratios for DO's setup are about 50% higher than for CDF's setup.) 
In both cases, the NLO corrections to the LO ratios are quite small, 12% or less. This 
suggests that the LO prediction of the ratio should also be fairly reliable for n = 4, though 
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confirmation awaits an NLO computation. In this case, however, the predicted ratio for 
n = 4 would lie significantly below the ratios for n = 2,3. We expect that non-perturbative 
effects will partly cancel in the ratios, making them theoretically more robust. Although 
the jet-production ratios can be used to give a rough estimate of higher-multiplicity jet total 
cross sections, for differential cross sections there is a strong dependence on the cuts and the 
number of jets. In particular, we found that the ratios depend strongly on the vector-boson 
Pt- These ratios of distributions have generic features which we described using a simple 
model capturing leading logarithms along with phase-space and parton-distribution-function 
suppression factors. The fits based on this model are surprisingly good and offer a simple 
parametrization of the theoretical predictions. 

Eventually, we would like to match the NLO results to parton showers and hadroniza- 
tion models, allowing non-perturbative effects to be modeled directly in an NLO program, 
instead of relying on LO-based tools to model them. This has been done for a variety of 
processes within the MC@NLO program [79| and the POWHEG method jsOj. It would 
be desirable to extend this matching to higher-multiplicity processes such as those presented 
in this present paper. A first step in this direction, linking BlackHat to an automated 
implementation of the FKS subtraction formalism used in MC@NLO, has been reported 
recently 



15 



81|. 



With the on-shell methods as implemented in BlackHat, we expect the computation of 
virtual corrections to cease presenting a bottleneck to obtaining new NLO results. A publicly 
available version of BlackHat is in preparation and is currently being tested in a variety 
of projects (see e.g. ref. Ml)- This version uses the proposed Binoth Les Houches interface 



for one-loop matrix elements [82[. It has been tested with both C++ and Fortran clients. 
We intend the public version to provide all processes that have been carefully tested with 
the full BlackHat code. The dipole-subtraction implementation that we used is available 



in the latest release of SHERPA 



NLO results will provide unprecedented precision for studies of vector bosons in associa- 
tion with many jets. They should prove useful for experimentally-driven determinations of 
backgrounds such as the invisible Z ^ vv background to missing-energy-plus-jets signatures. 
By measuring the corresponding Z — )■ Z"*"/" oxW^lv processes, with NLO predictions for 
cross-section ratios provi ding the necessary conversion factors, the Z vv background can 
be determined precisely [l9|, l20|. The results presented in this paper are examples of the 
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predictions that can be obtained using BlackHat in conjunction with SHERPA. We look 
forward to applying these tools to a wide range of studies of the forthcoming LHC data. 
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Appendix A: Kinematics and Observables 

In this appendix we give our definitions of standard kinematic variables used to charac- 
terize scattering events. The angular separation of two objects (partons, jets or leptons) is 

denoted by 



with the difference in the azimuthal angles, and Ay the difference in the rapidities. The 
rapidity is defined to be 



where E is the energy and pl is the component of the momentum along the beam axis (the 
z axis). The pseudorapidity 77 is given by 



DE-AC02-05CH11231. 




(Al) 




(A2) 




(A3) 



where 9 is the polar angle with respect to the beam axis. 
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Jets are formed using cluster or cone algorithms based on the angular separation (lAip . 
It is also possible to use the pseudorapidity in place of the rapidity. We have checked, for 
the production of Z, 7* + 3 jets at NLO, that using rj instead of y for the anti-/cr and kx 
algorithms with R = 0.5 makes no discernible difference on the cross section. 

The transverse energies of massless outgoing partons and leptons, Ex = ^pl + Py, can 
be summed to give the total partonic transverse energy, Ht, of the scattering process, 

Ht = Y.Et + Et + Ef . (A4) 
p 

All final-state partons p and leptons are included in Ht, whether or not they are inside jets 
that pass the cuts. As discussed in section UlIDl the variable Ht represents a good choice for 
the renormalization and factorization scale of a given event. Although the partonic version 
is not directly measurable, for practical purposes as a scale choice, it is essentially equivalent 
(and identical at LO) to the more usual jet-based total transverse energy, 

Ht = Y1 E't'j + E^^ + Ef . (A5) 
j 

The jet four-momenta are computed by summing the four-momenta of all partons that 
are clustered into them, 

pjr* = E • (A6) 

is jet 

The jet transverse momentum is then defined in the usual way, 

i^T = ^/iiW^Wy ■ (A7) 

Appendix B: Squared Matrix Elements at One Point in Phase Space 

As an aid to future implementation of Z, 7* + 3-jet production in other numerical codes, 
we present values of the one-loop virtual corrections to the squared matrix elements, day \ 
at one point in phase space. These contributions arise from the interference between the 
tree and one-loop amplitudes, summed over all colors and helicities, for Nc = 3 and Uf = 5 
massless quark flavors. As discussed in the text, we do not include the small effects from 
the top quark, or from axial or vectorial loop contributions (see fig. [2]). 

In table IIXI we present numerical values for three representative subprocesses. The other 
subprocesses are related to these by crossing symmetry or by change of coupling constants. 
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In the second line of table lIXt the presence of two identical quarks (after crossing all particles 
into the final state) implies that amplitudes are antisymmetrized under exchange of the two. 



day 




lA 




(ls2rf -> 3rf4^i5g6e-7e+) 


-8.3333333333 


-32.3745606495 


5.2374716277 




-8.3333333333 


-32.5180902998 


0.4387412994 


{lu'ig 3g4g5s6e-7e+) 


-11.6666666667 


-42.3279266518 


-15.2326082853 



TABLE IX: Numerical values of the normalized virtual correction to the squared matrix elements, 
day\ at the phase-space point given in eq. (|B3p . for three basic partonic subprocesses for Z+ 3-jet 
production at a hadron collider. We give both the finite parts and the coefficients of the poles in e. 

We quote numerical results for the ultraviolet-renormalized virtual corrections in the 
't Hooft-Veltman variant of dimensional regularization 8^. The remaining singularities in 
the dimensional regularization parameter e = {4 — D)/2 arise from the virtual soft and 
collinear singularities in the one-loop amplitudes. 

The quoted values are for the ratio of the virtual corrections to the tree-level squared 
matrix element da^^\ Explicitly, we define the ratio. 



da 



(1) 



da 



(1) 



(Bl) 



^ 87ra5Cr(e)da(o) ' 

where we have also separated out the dependence on the strong coupling as and the overall 
factor cr(e), defined by 



(47r)2-^ r(l - 2e) 

The coupling constants, mass and width of the Z boson are given in section IIII A[ 
We choose the phase-space point given in eqs. (9.3) and (9.4) of ref. |25| . 



h 



2 (1, — sin 9, — cos 9 sin 0, — cos 9 cos ( 
^ (1, sin 9, cos 9 sin 0, cos 9 cos 0) , 
|(1, 1,0,0), 
— (1, cos/3, sin/3, 0) , 



10 
12 

ki + k2- ks 



1, cos a cos /3, cos a sin 13, sin a) 
1, cos 7 cos /3, cos 7 sin /3, sin 7) , 
ki — k^ — kg , 



(B3) 
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where 

^=4' ^^=6' "=3' ^ = T' ^°^^ = -128' ^^'^ 

and the renormahzation scale fiji is set to fiR = 7 GeV. We have flipped the signs of ki and 
^2 compared to ref. |25], to correspond to 2 — )■ 5 kinematics, instead of — ?■ 7 kinematics. 
The labehng of the parton and lepton momenta is indicated exphcitly in the first column of 
Table IXl 
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